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The structure of the second 2 + resonance in 6 Li is investigated with special emphasis on its 
isospin components. The wave functions are computed in a three-body model (a+n+p) using 
the hyperspherical adiabatic expansion method combined with complex scaling. In the decay into 
three free particles the symmetry conserving short-range interaction dominates at short distance 
whereas the symmetry breaking Coulomb interaction dominates at intermediate and large distances 
resulting in substantial isospin mixing. We predict the mixing and the energy distributions of the 
fragments after decay. Computations are consistent with available experiments. We conjecture that 
nuclear three-body decays frequently produce such large isospin mixing at large distance where the 
energy distributions are determined. 
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Introduction. The spatial extension of halo states de- 
pends sensitively on their binding energy |lj . The reason 
is that the outer particles forming the loosely bound halo 
are confined by an attraction or a barrier of moderate 
size. This sensitivity to the energy must be particularly 
important for resonances owing their existence to a con- 
fining barrier. For the same reason, due to the effect of 
the Coulomb interaction, the structure of isobaric ana- 
log resonances can be very different. This influence can 
also be reflected in mixing of different isospin compo- 
nents. For ordinary states the isospin mixing has been 
established to be in the range 10 _4 -1(T 5 0. For isobaric 
analogs to halo states the mixing could be substantially 
larger 

Early theoretical investigations of decay of three-body 
analog halo states in light nuclei indicated that isospin is 
very well conserved 0, Q . A very small isospin mixing is 
also deduced from the reduced measured branching ratio 
of decay of the highest known 2 + resonance in 6 Li (with 
excitation energy of 5.37 MeV and associated to isospin 
1) into the a-deuteron two-body system [f| 0] (admix- 
ture of 8 • 10 -3 or less). Nevertheless, isospin mixing in 
6 Li-states is also investigated in a number of reaction ex- 
periments 2 H(cy, a) 2 H* where 2 H* is the neutron-proton 
system in a spin singlet relative state which cannot be 
populated when isospin is conserved. Strong evidence 
is claimed for the formation of 2 H* which means sub- 
stantial isospin mixing and especially when the incident 
energy is low [1, [t| [Ic| • Reasonable agreement between 
measurements and simple model calculations can only be 
obtained with inclusion of both spin zero neutron-proton 
s-wave ( 2 H*) and nucleon-a relative d- waves. Some of 
these authors § suggest the large isospin mixture to 
be due to direct reactions bypassing the 6 Li resonance, 
whereas others @ call for better methods for including 
the Coulomb interaction in three-body calculations. 

Recent investigations emphasized that the structure 
of three-body resonances can change substantially from 



small to large distances [ll|, [l2|, [l3|. This means for 
instance that the amount of isospin may vary substan- 
tially with the (hyper) radial coordinate describing the 
wave function. In other words, the relative weights of 
different partial wave components can be very different 
from small to large distances, e.g. in 6 Li (a + n + p) 
for neutron and proton in relative s-waves and in triplet 
or singlet spin states corresponding to T = or T = 1, 
respectively. The absence of decay of the 2 + resonance 
in 6 Li into the deuteron channel is not a direct evidence 
for a correspondingly small admixture of isospin 0. The 
dominant isospin 1 wave function at small distance can 
decay into continuum final states with isospin provided 
the coupling is sufficiently strong. This is most likely at 
distances just outside the ranges of the short-range inter- 
actions where the Coulomb interaction still is substantial 
and completely dominating. 

The purpose of the present letter is to investigate the 
isospin conservation in the 2 + three-body resonance in 
6 Li, that is the isobaric analog state of the correspond- 
ing 2 + resonance in 6 He and 6 Be, for which isospin zero 
components are forbidden. We start with a sketch of the 
necessary theoretical framework. We then discuss the 
structure of the resonance and present the final state en- 
ergy distributions. We finish with a brief summary and 
the conclusions. 

Theoretical framework. Bound states and resonances 
are obtained with the hyperspherical adiabatic expansion 
method combined with complex scaling. These three- 
body wave functions \& are then expressed as a linear 
combination of the complete set of functions {<& n (p, f2)} 
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where p is the hyperradius and the five hyperangles fl = 
{a, H, x ,H,y} can formally be chosen in any of the three 
Jacobi sets, or as we prefer the wave function can be given 
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in terms of three Faddeev components each associated 
with one Jacobi set. Each of these Faddeev components 
are in turn expanded in partial waves with a basis of 
corresponding hyperspherical harmonic functions. 

The functions <£>„ (p, ft) are the eigenfunctions of the 
angular part of the Faddeev equations, and the radial 
coefficients f n (p) are obtained from the coupled set of 
radial equations where the eigenvalues of the angular part 
enter as effective adiabatic potentials [T3 |. Due to the 
complex scaling the wave functions fall off exponentially 
at large hyperradii for bound states and resonances with 
width-to-energy ratios less than twice the rotation angle. 

The components of the solution directly contain infor- 
mation about the symmetry and hence about the isospin 
mixing, which can vary with hyperradius. This point 
has apparently not been appreciated in the old analyses 
0, 0). For instance, in 6 Li (n+p+a) the total isospin 
is obtained by coupling neutron and proton isospins to 
either or 1. Components with both isospins are in 
principle permitted, like for instance (£ x ,£ y , L, S, J) — 
(0, 2, 2, 0, 2), (1, 1, 2, 0, 2), where £ x and £ y are orbital an- 
gular momenta related to the Jacobi coordinates x and 
y (with x proportional to the neutron-proton distance), 
and L, S and J are the total orbital, spin and total an- 
gular momenta. In fact these two components do con- 
tribute, and are coupled due to the presence of the a- 
particle, which Coulomb interacts with the proton but 
not with the neutron. This can totally break the (isospin) 
symmetry if the symmetry breaking interaction is dom- 
inating as it is for large distances. Thus, small isospin 
mixing at small distance can be compatible with a very 
large isospin mixing at large distance. 

The kinetic energy distribution of the fragments after 
decay of a resonance is, except for a phase-space factor, 
obtained as the absolute square of the total wave function 
in coordinate space for a large value of the hyperradius 
p, but where the five hyperangles are interpreted as in 
momentum space (l2l fl3| . After integration over the four 
hyperangles (tt x ,£l y ) describing the directions of the two 
Jacobi momenta, k x and fc y , conjugate to x and y, the 
probability distribution as function of k y oc cos 2 a, where 
a is the fifth momentum hyperangle, is given by 

P(fc 2 ) ocP(cos 2 a) ocsin(2a) J dQ x dSl y \^(p, a,n x ,Q y )\ 2 . 

The kinetic energy of the third particle is proportional to 
k y oc cos 2 a which then gives the energy of the particle 
relative to its maximum possible energy in the decay pro- 
cess. These observables carry information about initial 
state and decay mechanisms. 

Details of the calculations. The 2 + resonance in 6 Li 
has been computed using the same a-nucleon interac- 
tion as in [la ] for 6 He. For the neutron-proton poten- 
tial we use the one in Hi'. Components with relative 
two-body orbital angular momenta up to 4 are consid- 
ered. The main components included in the calculation 
are shown in table U (the first column labels the compo- 
nents). A proper choice of the maximum value of the 



TABLE I: Components included for the 2 + -state in 6 Li. The 
left part refers to the first Jacobi set (x from neutron to pro- 
ton), and the right part to the second and third Jacobi sets (x 
from one of the nucleons to the a-particle). The first column 
numbers the different components. 
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FIG. 1: Outer panel: Real parts of the (complex scaled) ef- 
fective adiabatic potentials for the 2 + resonance in 6 Li. Inner 
panel: Real parts of the radial wave functions associated to 
the three most contributing effective potentials (indicated by 
the thick curves in the outer part). 



hypermomentum (K max ) for each of them is crucial to 
obtain a correct convergence of the effective potentials 
at a sufficiently large distance. The K max value for the 
components 6 to 11 in the left part of table [J is relatively 
large to ensure an accurate calculation of their contribu- 
tion, since these are precisely the components with zero 
isospin (T in the table) that are not allowed in 6 He or 
6 Be. For the remaining components (not shown in the 
table) the -Km^-value is at least 20. 

The outer part of FigQ]shows the real parts of the effec- 
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tive potentials obtained after use of the (complex scaled) 
hyperspheric adiabatic expansion method. The lowest ef- 
fective potential converges towards the deuteron binding- 
energy, and appears due to the inclusion of the compo- 
nents with zero isospin in the neutron-proton channel. 
These potentials are indistinguishable from the ones ob- 
tained when the basis size is reduced by a factor of two. 
This fact guaranties the convergence of the potentials at 
least up to 100 fm. 

Calculation of the resonance wave function requires 
specification of the corresponding boundary condition. 
As shown in T7j], a simple box boundary condition 
at a sufficiently large distance is enough to obtain a 
resonance wave function with the proper asymptotics. 
This requires extrapolation of the effective potentials for 
p values beyond 100 fm. An expansion of the form 
A/ p + B I p 2 + C I p 3 + ■ ■ ■ is used (except for the low- 
est one). The non-adiabatic coupling functions P n n'{p) 
and Qnn'(p) (see [l4|) are extrapolated as Q n n{p) = 
A Q n Jp 2 + B Q,m/p 3 + c Q n n/p 4 f° r diagonal Q's (diago- 
nal P's are zero), and as Apg/p 3 + Bp^/p 4 + Cp^/p 5 
for the non-diagonal P's and Q's. 

A box boundary condition at p mQa: =1000 fm gives 
rise to a 2 + resonance in 6 Li with energy and width 
(E R ,T R )=(1.67,0.51) MeV (5.37 MeV excitation en- 
ergy), that agrees with the experimental value of 
( J B J? ,r i? ) = (1.67±0.02,0.54±0.02) MeV [3 (the reso- 
nance energy is indicated in FigQ] by the dashed line). 
However, the extrapolation used for the effective poten- 
tials implies that the asymptotics of the radial wave 
functions must be F^{rj,Kp) — iG^{rj, up), where n = 
y / 2mE/h 2 , and are the regular and irregular 
Coulomb functions, and the Coulomb charge r\ and the 
index £ can be easily obtained from A, B, and -Ag nn . 
When this asymptotic condition is imposed, a much 
smaller value of p m ax is enough to obtain the resonance 
wave function (but still p max >100 fm, and the expan- 
sions of the effective potentials, P's, and Q's are re- 
quired). 

In the inner part of FigQ] we show the real parts of 
the computed radial wave functions. A complex scaling 
angle of 0.10 rads has been used. To make the picture 
cleaner we only show the ones associated with the three 
most contributing effective potentials (thick curves in the 
outer part of the figure). It is important to note that the 
contribution from the lowest adiabatic potential, the one 
holding deuteron at large distances, is very small, provid- 
ingabout 10 _3 % of the norm, which is consistent with 
[sul- It is also necessary to investigate whether the ra- 
dial wave functions have already reached the asymptotic 
behaviour for p < 100 fm, i.e., in the p-region where 
the calculation is purely numerical. This would mean 
the asymptotics given by the extrapolated effective po- 
tentials, P's, and Q's, is consistent with the numerical 
results. 

This is tested in Fig [21 where the real and imagi- 
nary parts of the three most relevant radial wave func- 
tions (thick curves) are shown and compared to the ex- 




FIG. 2: Thick curves: Radial wave functions corresponding to 
the three most contributing effective potentials. Thin curves: 
Expected asymptotics (F^(ri,Kp) — iG^(ri, up)) from the ex- 
trapolated potentials for p >100 fm. 



pected asymptotics as given by the regular and irregular 
Coulomb functions, with Coulomb charges and indices 
computed numerically from the extrapolation of the po- 
tentials (thin curves). The matching between the numer- 
ical wave functions and the asymptotics is already very 
good at about 60 fm, clearly below the p-value (100 fm) 
from which the extrapolations are used. The labels 1, 2, 
and 3 refer to the deepest thick potential, intermediate 
thick potential, and repulsive thick potential in the outer 
part of FigfTJ 

From the three-body wave function in Eq.([T]), we define 
the total weight (as function of p) as: 

W(p) — J sin 2 acos 2 adadfl x dfl y \^(p, a, fl x , Q y )\ 2 , 

(3) 

When writing VP in the first Jacobi set [x from neutron 
to proton) we find, after integration over p in Eq.([3]), 
that roughly 82% of the weight is given by the three first 
components in the left part of table Q] (38%, 20%, and 
24%), while the remaining 18% is distributed among the 
other components. 

It is remarkable that the last components from 6 to 11 
in the left part of table U (associated to zero isospin) ac- 
cumulate about 4.5% of the integrated weight, most of it 
corresponding to component 9. This is due to the third 
adiabatic potential (repulsive thick potential in Fig[T|), 
whose corresponding eigenfunction $3(p, f2) is dominated 
at intermediate distances by the components with zero 
isospin in the neutron-proton channel. This is shown 
in the upper part of Fig[31 where we show, as a func- 
tion of p, the main contributions to $3(p, f2) from the 
components in the left part of table HI The thick curves 
correspond to components 8 and 9 in the table (T=0). 
As observed in the figure, these components have a non- 
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FIG. 3: Top: Main contributions (as a function of p) of the 
components in the left part of table [T] to the eigenfunction 
$3 (p, fi) associated to the third adiabatic potential (thick re- 
pulsive potential in Fig[T}. The thin and thick curves cor- 
respond to components associated to T=l and T—0, respec- 
tively. Bottom: Weight of the T=l (dashed) and T=0 (dot- 
dashed) components in the 2 + resonance wave function in 6 Li 
(scaled by a factor 2.5). The thick and thin solid curves are 
the relative weight of the T—0 components and the relative 
weight of the 9 th component in the left part of table HI respec- 
tively. 



negligible weight at intermediate distances. In particu- 
lar, component 9 gives a large contribution from 20 to 
50 fm. Beyond 50 fm component 3 dominates (T=l), 
but still a 5% contribution from component 9 is present. 
The rapid transition in $3 from component 9 to 3 reflects 
that the isospin symmetry is totally broken, because only 
the Coulomb interaction is active. The lowest centrifugal 
barrier with £ y =l is then abruptly preferred over t y =2. 
At short distances, as seen on the bottom of Fig(2l the 
radial coefficient corresponding to the third adiabatic po- 
tential is negligible, but beyond 20 fm, the amplitude of 
/3 (p) is similar to the one in /1 (p) and fa (p) . 

Let us denote now by Wt=i and Wt=o the contribu- 
tions to W(p) in Eq.Q from the T=l and T=0 com- 
ponents in the three-body wave function, respectively. 
These contributions are shown (scaled by a factor 2.5) 
by the dashed and dot-dashed curves in the bottom part 
of Figj3j At short distances the T=l contribution clearly 
dominates, while Wt=o becomes relevant at intermediate 
p's. This is more clearly seen by the thick solid curve, 
that shows the relative weight of the T=0 components. 
From 20 to 50 fm this weight reaches up to 40% of the 
total. This region coincides with the one where $3 (p, Q) 
has a relevant contribution from the T=0 components 
(upper part of the figure). Beyond 50 fm Wt=o stabi- 
lizes at about 10%. In the figure the thin curve shows the 
relative contribution to the total weight from component 
9. This contribution gives most of the T—0 contribution, 
and governs the general behaviour of Wt=q/W. 
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FIG. 4: Thick curves: Energy distribution of the a (solid), 
neutron (dashed), and proton (dot-dashed) after decay of the 
2 + resonance in 6 Li. Thin curves: The same energy distribu- 
tions after excluding the T—0 components. 



The non-negligible isospin mixing found at large dis- 
tances can be relevant for those observables sensitive to 
the asymptotic behaviour of the wave function, like for 
instance the energy distributions of the fragments after 
decay. In Fig 0] we show the a (solid), neutron (dashed), 
and proton (dot-dashed) energy distributions according 
to Eq.([2]). The thick curves show the total distributions, 
while the thin ones give the same distributions when the 
T=0 components have been excluded. The results shown 
in the figure are very stable for p values in Eq.([2]) rang- 
ing from 65 to 85 fm. In particular the curves shown in 
the figure correspond to p — 75 fm. As seen in the fig- 
ure, inclusion of the T=0 components produce a visible 
change in the energy distributions, although the general 
behaviour of the distributions does not change. 

Summary and conclusions. We have investigated dy- 
namic isospin mixing in nuclear resonances. We illustrate 
by the highest known 2 + resonance of the three-body 
system 6 Li, where components with different isospin si- 
multaneously can be present. The Coulomb interaction 
between a-particle and proton is breaking the isospin 
symmetry and mixing isospins of and 1. The isospin 
zero components essentially only appear in the neutron- 
proton continuum. The deuteron is populated in the de- 
cay by about 10~ 3 %, that is consistent with the exper- 
imental upper limits [1, 0] • The amount of isospin zero 
is consistent with that found in experiments [3, Q (up 
to 30% isospin mixing for low energies), while here the 
isospin mixing is due to the decay process and not direct 
reactions as suggested by these authors. 

We have used the complex scaled, hyperspheric adia- 
batic expansion method with an extraordinary large basis 
for the components with zero isospin. The isospin con- 
tent of the accurately computed resonance wave functions 
vary substantially from small to large distances. The rel- 
ative isospin contribution is small at small distances 
where the main contribution resides. This relative con- 
tribution reaches about 40% at intermediate distances, 
and stabilizes beyond 50 fm at roughly 10% of the to- 
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tal. The total contribution integrated over all distances of dynamic isospin mixing is a common feature in decays 
of the isospin components is about 4% which is much of nuclear three- (or more-) body resonances, 
larger than for ordinary stable nuclei. This mechanism 
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